Local density of states of a strongly type-II d-wave superconductor: 
The binary alloy model in a magnetic field 
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We calculate self-consistently the local density of states (LDOS) of a d-wave superconductor 
considering the scattering of the quasiparticles off randomly distributed impurities and off externally 
induced vortices. The impurities and the vortices are randomly distributed but the vortices are 
preferably located near the impurities. The increase of either the impurity repulsive potential or 
the impurity density only affects the density of states (DOS) slightly. The dominant effect is due 
to the vortex scattering. The results for the LDOS agree qualitatively with experimental results 
considering that most vortices are pinned at the impurities. 
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Experimental evidence suggests that the pairing sym- 
metry in high temperature superconductors (HTSC) is 
d-wavei A good description of the non-conventional su- 
perconducting phase is obtained using a standard BCS 
approach but a clear understanding of the normal phase 
of these materials remains a hard challenge. If we ap- 
ply a strong enough magnetic field, these materials (be- 
ing strongly type-II superconductors) will enter a vortex 
phase. Scanning Tunneling Microscopy (STM) studies of 
HTSC have revealed a very different quasiparticle struc- 
ture from that predicted by the pure d-wave BCS single 
vortex models*. In a pure d-wave there are extended gap- 
less states, a fourfold symmetric shape of the LDOS and a 
zero-bias conductance peak. Experimentally, however, it 
is observed: i) absence of zero energy peaks, ii) absence of 
coherence peaks close to the vortex, iii) low energy states 
with an energy ~ 5.5meV for YBCO and ~ 7.7meV for 
BSCCO, and iv) absence of 4- fold symmetric star-shaped 
LDOS*. The coherence peaks are recovered about lfjA 
from the core center— and the core states are localized 
decaying exponentially with distance (~ 22 A}-. Possi- 
ble reasons for the failure of the pure d-wave theory to 
explain the experimental results have included a mixed 
pairing of the type d x 2_ y 2 + id xy , considering that the 
vortices have antiferromagnetic cores such that localized 
magnetic order coexists with superconductivity or charge 
order fluctuations*!*. 

Recently, a pure d-wave pairing for a vortex lattice has 
yielded results that are in good qualitative agreement 
with experiments*. Indeed a calculation of the local den- 
sity of states shows that close to the vortex position the 
coherence peaks disappear, there are significant low en- 
ergy peaks and there is no zero energy enhancement of 
the density of states*. However, in most systems disor- 
der is present for instance in the form of impurities. The 
presence of the impurities affects the quasiparticles (QP) 
in two ways: the QP scatter off the impurities due to 
potential scattering (if they are nonmagnetic) and the 
impurities pin the vortices also affecting the density of 
states of the QP, particularly at low energies. 

The separate effects of the scattering of the quasiparti- 



cles from the impurities and from the vortices have been 
studied before. The details of the impurity disorder are 
relevant and a consistent picture of the various possible 
scenarios has been obtained*. Studies of the superfluid 
stiffness due to the presence of the impurities have re- 
vealed that, even though it gets lower, the decrease is 
smaller than expected. The reason is that the order pa- 
rameter is only significantly affected very close to the 
impurities and largely unaffected elsewhere. Therefore 
the order parameter is very non-homogeneous and a fully 
self-consistent calculation is required*. Considering ex- 
clusively the effect of the scattering off homogeneously 
distributed vortices (no impurities) it has been shown 
that in the lattice case the low energy states are ex- 
tended Bloch states^ instead of Dirac Landau levels**. 
Also, it was shown that the quasiparticles besides feeling 
a Doppler shift caused by the moving supercurrentsi*, 
also feel a quantum "Berry" like term due to a half-flux, 
4>o/2, Aharonov-Bohm scattering of the quasiparticles by 
the vortices. The effect of a random vortex distribution 
was considered taking random and statistically indepen- 
dent scalar and vector potentials**. A finite density of 
states was predicted at zero energy. Also, considering 
randomly pinned discrete vortices the density of states 
was calculated displaying low energy peaks and no co- 
herence peaks** but at low energies a power law devi- 
ation from a finite zero energy value was found, where 
both the zero energy value and the exponent depend on 
the magnetic field and on the Dirac anisotropji*. Also, 
it was found that even though the low energy states are 
strongly peaked at the vortex cores, they appear to re- 
main extended. An approximate scaling of the density of 
states was found at low energies*. 

In this work we study the combined effects on the quasi- 
particle spectrum of the scattering off impurities and vor- 
tices induced by an external magnetic field. We model 
the disorder using the binary alloy model** where the im- 
purities are distributed randomly over the lattice sites. 
At each impurity site it costs an energy U to place an 
electron (it acts as a local shift on the chemical poten- 
tial). The impurities are randomly distributed over a 
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L x L periodic two-dimensional lattice and play the role 
of pinning centers for the vortices. It is favorable that a 
vortex is located in the vicinity of an impurity 15 . Taking 
into account a given distribution of the positions of the 
impurities {rf N , the distribution of the positions of 
the vortices {r"} i=1 N is chosen in such a way it mini- 
mizes the total vortex energy given by £ = £ v +£ v , where 
£ v = U v J2 r * r » Kq (\r^ — r v : I /A) is the repulsive interac- 
tion energy between the vortices in the London regime 
and £ p = K p J2 r ^,r p f (| r i _ r j I / r p) is tne P mn ing en- 
ergy associated with the impurities acting as pinning cen- 
ters for the vortices. The interaction between the vortices 
is not significantly screened since the penetration length 
is very large. In the equations above U v — (0o/47tA) 2 
is the energy of interaction between two vortices, Kq(t) 
stands for the zero order Hankel function, U p < is the 
pinning strength created by an impurity and f(r/r p ) is a 
rapidly decreasing function for r/r p > 1. In our model we 
assume that the pinning energy is much larger than the 
interaction between vortices \U P | 3> IA V and r p ~ 8 where 
5 is the lattice constant. In that case, as N v <C N p each 
vortex is preferably pinned in the close vicinity of an im- 
purity. As we take the London limit, which is valid for low 
magnetic field and over most of the H — T phase diagram 
in extreme type-II superconductors such as cuprates, we 
assume that the size of the vortex core is negligible and 
place each vortex core at the center of a plaquette. So in 
the limit of strong pinning described above, each vortex 
will be pinned in the center of one of the four plaquettes 
surrounding a site hosting an impurity. The plaquetes 
selected by the vortices are those minimizing the inter- 
action energy £ v between the vortices. 

Once the impurity positions {r^} i=1 N are fixed and 
the correlated vortex positions {r" } i=1 N are determined 
we are able to calculate de quasiparticle spectrum us- 
ing the BdG equations H(r)^> n (r) = e n ^ n (r), where 
^(r) = (m* (r), v* (r)). It is convenient* 4 ^ to perform 
a unitary gauge transformation. After carying out this 
gauge transformation the Hamiltonian reads 

-1) « 

where 

s 

A = J2 e lAs{r)+nS «A(r, r+S)s s . 
s 

The phase factors are given by Vg (r) = f* +S k^ • dl 

and As{r) = \ f* +S (k^(l) - kf (1)) • efl where the vector 
frk^(r) = mvj(r) = ftV^ — (e/c)A(r) is the superfluid 
momentum vector of the effective /i = A, _B-supercurrent. 
This quantity can be calculated for an arbitrary configu- 



ration of vortices like 

%w- 27r y (27r )2 k 2 + A - 2 2^ e ■ w 

Here A is the magnetic penetration length and the sum 
extends over all vortex positions. A(r) is the vector po- 
tential associated with the uniform external magnetic 
field B = V x A. The vortices A are only visible to 
the electrons and the vortices B are only visible to the 
holes. Each resulting /x-subsystem is then in an effective 
magnetic field B^ ff = -Sf V x = B-0 o z£ 4 6(r-r? ) 
where each vortex carries now an effective quantum mag- 
netic flux </>o . For the case of a regular vortex lattic e) 4 ! 16 , 
these effective magnetic fields vanish simultaneously on 
average over a unit magnetic cell containing two vortices, 
one of each type. More generally, in the absence of spa- 
tial symmetries, as it is the case for disordered systems, 
these effective magnetic fields B^ j4 ' s vanish if the num- 
ber of vortices of the two types are in equal number, 
./V^ = Nb , and equal to the number of elementary quan- 
tum fluxes 0o °f the external magnetic field penetrat- 
ing the system. The sums are over nearest neighbors 
(S = ±x, ±y on the square lattice) and the operator sg 
is defined through its action on space dependent func- 
tions, §su(r) — u(r + S). The energy U{r) is the impu- 
rity potential which takes the value U > at the sites 
{ r f n hosting an impurity and zero elsewhere. The 
operator f]$ = e™ Sy ss carries the symmetry of the d- 
wave order parameter. The disorder potential induced 
by the impurities and the inhomogeneous superfluid ve- 
locities induced by the vortices strongly affect the pairing 
potential A(r,r+£) defined over the link [r, r+S] of the 
two-dimensional lattice. Thus for a given configuration 
of the impurity positions and of the vortex positions the 
pairing potential A(r, r+<5) is calculated self-consistently 
until convergence is obtained for each individual link. On 
each lattice site we can define the amplitude of the rf-wave 
order parameter as A4 (r) = l)" 5 " A(r, r+^). The 

amplitude of the order parameter is strongly suppressed 
in the vicinity of impurities and vortices. 

As the effective magnetic fields experienced by the par- 
ticles and the holes vanish on average, within the gauge 
transformation we are allowed to use periodic boundary 
conditions on the square lattice (ty(x + nL,y + mL) = 
*&(x,y) with n,m £ Z). The L x L original lattice be- 
comes then a magnetic supercell where the impurities are 
placed at random and where the vortices are placed in 
such a way as to minimize their total energy. The dis- 
order induced by the impurities in the system is then 
established over a length L. Thus in order to com- 
pute the eigenvalues and eigenvectors of the Hamilto- 
nian Q we seek for eigensolutions in the Bloch form 
^+ k (r) = e- lkr (C/* k ,y„* k ) where k is a point of the 
Brillouin zone. We diagonalize then the Hamiltonian 
e" We for a large number of points k in the Bril- 
louin zone and for many different realizations (around 
100) of the random impurity positions and of the corre- 
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FIG. 1: Quasiparticle density of states for a 20 x 20 sites 
system with 4% of impurities, with U = 5t, £f = 1.2 and 
V = — 2.3i. We consider three cases: no vortices (T), 4 vor- 
tices pinned at locations ensuring the minimization of the 
total vortex energy (□), and 4 vortices located at random 
(•). For each case the density of states is averaged over 100 
configurations. 




FIG. 2: Quasiparticle density of states for a 20 x 20 sites. 
Left panel: system with 2% of impurities, 4 vortices pinned at 
locations ensuring the minimization of the total vortex energy 
and U = 2*(o), 5*(B), 10t(D) and 100t(^). The case of 4 vor- 
tices pinned at random in a system without impurity is also 
presented (•). Right panel: system with 4 vortices pinned at 
locations ensuring the minimization of the total vortex energy, 
U — 5t and the following percentages of impurity: 0% [vor- 
tices pinned at random] (•), l%(o), 2%(B), 4%(D) and 6%0). 



lated vortex positions. 

In Fig. ^we show the density of states averaged over 
100 configurations for a moderate impurity concentra- 
tion and for a typical value of U = 5t. We compare 
the cases with no vortices and with a low vortex density 
(considering both an energetically favorable distribution 
of the vortices and a fully random configuration). The 
self-consistent calculation of the order parameter gives a 
maximum amplitude, Ao, of the order of 0.5i. It is clear 
that when there are no vortices present, coherence peaks 
appear at e ~ Aq. If the vortices penetrate the sample 
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FIG. 3: Quasiparticle local density of states for different 
points on the lattice. The impurity concentration is 0.75%. 
We have chosen a particular distribution of the impurities and 
vortices where we can find the various possibilities of a site 
far from any impurity or vortex (Bulk), an impurity site with 
no vortex nearby (Imp), a site close to a vortex and no impu- 
rity nearby (Vx) and an impurity site with a vortex attached 
(Imp+Vx). For each of these 4 possibilities we have calcu- 
lated the quasiparticle local density of states along a path of 
six numbered sites 1, . . . , 6. On the top most pannel the 4 
particular paths are shown on the top of the corresponding 
profile of the d-wave order parameter Ad(r). The path located 
on the center of the lattice corresponds to the Bulk case, the 
one centered on the site (6,12) corresponds to the Imp+Vx 
case, the one centered on the site (3,3) corresponds to the Vx 
case, and the one centered on the site (11,3) correspond to 
the Imp case. 
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these peaks disappear. Without vortices the density of 
states vanishes at zero energy as found before and if we 
include the magnetic field the density of states becomes 
finite 6 . Also it is clear that if the vortices are fully ran- 
domly distributed the DOS is larger than the one where 
the vortices tend to be pinned at the impurity sites. The 
results are therefore qualitatively similar to the ones ob- 
tained when there are no impurities present (U = 0) 6 . 
In Fig. |2 we show the influence of the impurity concen- 
tration and of the repulsive local potential, U. Changing 
the impurity concentration leads to no qualitative differ- 
ence except that there is a slight increase in the DOS. 
The same happens with the change of U. Both sets of 
results indicate that the strong effect is the scattering off 
the vortices. 

More detailed information about the scattering of the 
quasiparticles is obtained calculating the LDOS for a 
given impurity /vortex configuration. It is defined by 

p(r, e) = J2 (KMI 2 + K(r)\ 2 ) 6(e - e n ). 

n 

In Fig. |31 we compare the LDOS at four different sites, 
i) At a site in the bulk the band-structure profiles are 
somewhat similar to the case of a vortex lattice^. The 
coherence peaks are evident, the zero energy density of 
states is very small (but not strictly zero), and the low 
energy peaks are smeared out. ii) At an impurity site 
(and no vortex nearby) the same structure is apparent 
except that since the impurity potential is repulsive the 
density of states at the impurity site is considerably de- 
pleted at low energy (for instance in the very large U 
limit the density of states is virtually zero at the impu- 
rity site), iii) In the vicinity of a vortex site (but far 
from any impurity) the structure is very similar to the 
case obtained before^ with no coherence peaks close to 
the vortex and an enhanced zero energy density of states 
near the vortex, iv) Finally the interesting case of a lo- 



cation where a vortex is bound to an impurity reveals 
that the coherence peaks are recovered very close to the 
vortex. Also, the low energy density of states at the im- 
purity site is increased with respect to the (Imp) case, 
due to the vortex nearby. However the density of states 
is considerably smaller than for the case (Vx) of a vortex 
far from any impurity. 

The results obtained previously for the vortex lattice 
case with no impurities explain qualitatively the DOS re- 
sults but are not realistic. In the experimental systems 
disorder is present and its effect must be taken into ac- 
count. The results show that the dominant effect on the 
quasiparticle DOS is due to the vortex scattering. The 
presence of an impurity basically renormalizes the DOS 
except when the impurity is strongly repulsive where the 
density of states is significantly depressed near the impu- 
rity. This is the unitary limit where a gapped system is 
predicted in the absence of magnetic fieldi The quan- 
tum effect originated in the Aharonov-Bohm scattering 
of the quasiparticles circulating around a vortex line has 
been shown to have considerable effects on the density of 
statesiSiii. Significant changes with respect to the clas- 
sical Doppler shift effectii occur at low energies^. The 
results obtained in this paper show that the addition of 
impurity scattering is not very significant and the Berry 
phase is dominant, as argued before^. The results are 
very similar to the ones obtained experimentally with 
STM in ref£, except for the finite density of states at 
zero energy in the vicinity of an isolated vortex. The ex- 
perimental results are therefore more consistent with the 
situation where all or most of the vortices are pinned to 
the impurity sites. At these points, even though there is 
an enhancement of the low energy density of states with 
respect to an impurity with no vortex attached, the in- 
crease is reduced by the presence of the impurity with 
respect to the case of a vortex but no impurity nearby. 
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